Equation of State for Strange Quark Matter in a Separable Model 



C. Go eke 

Fachbereich Physik, Universitdt Rostock, D-18051 Rostock, Germany 



D. Blaschke 

Institute for Nuclear Theory, University of Washington, Box 351550, Seattle, WA 98195 
Fachbereich Physik, Universitdt Rostock, D-18051 Rostock, Germany 
Bogoliubov Laboratory for Theoretical Physics, JINR, RU-14-1980 Dubna, Russia 

o 

A. Khalatyan 

Department of Physics, Yerevan State University, 375025 Yerevan, Armenia 



, H. Grigorian 

Fachbereich Physik, Universitdt Rostock, D-18051 Rostock, Germany 
Department of Physics, Yerevan State University, 375025 Yerevan, Armenia 



CO 
00 



o 



X 

5-H 



Abstract 



We present the thermodynamics of a nonlocal chiral quark model with sep- 
Qh. arable 4-fermion interaction for the case of U(3) flavor symmetry within a 

r- 1 \ functional integral approach. The four free parameters of the model are fixed 

by the chiral condensate, and by the pseudoscalar meson properties (pion 
mass, kaon mass, pion decay constant). We discuss the T = equation 
of state (EoS) which describes quark confinement (zero quark matter pres- 
sure) below the critical chemical potential [i c = 333 MeV. The new result of 
the present approach is that the strange quark deconfmement is separated 
from the light quark one and occurs only at a higher chemical potential of 
fj, CtS = 492 MeV. We compare the resulting EoS to bag model ones for two 
and three quark flavors, which have the phase transition to the vacuum with 
zero pressure also at fi c . 

We study quark matter stars in general relativity theory assuming (5- 
equilibrium with electrons and show that for configurations with masses close 
to the maximum of stability at M = 1.62 -j- 1.64 M & strange quark matter 
can occur. 

PACS numbers: 04.40.Dg, 12.38.Mh, 26.60. +c 
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I. INTRODUCTION 



Since the discovery of the parton substructure of nucleons and its interpretation within 
the constituent quark model, much effort has been spent to explain the properties of these 
particles. The phenomenon of confinement, i.e. the property of quarks to exist only in bound 
states as mesons and baryons in all known systems, poses great difficulties for a describing 
theory. So far the problem has been solved by introducing a color interaction that binds all 
colored particles to "colorless" states. 

However, it is believed and new experimental results underline it, that at very high 
temperatures exceeding 150 MeV, or densities higher than three times nuclear matter den- 
sity, a transition to deconfined quark matter can occur. Besides of the heavy ion collisions 
performed in particle physics, the existence of a deconfinement phase and its properties is of 
high importance for the understanding of compact stars @ in astrophysics. These, that are 
popular as Neutron stars, imply core densities above three times the nuclear saturation den- 
sity so that quark matter is expected to occur in their interior || and several suggestions 
have been made in order to detect signals of the deconfinement transition ||||. 

Unfortunately, rigorous solutions of the fundamental theory of color interactions (Quan- 
tumchromodynamics (QCD)) for the EoS at finite baryon density could not be obtained yet, 
even Lattice gauge theory simulations have serious problems in this domain [^]. To describe 
interacting quark matter it is therefore necessary to find approximating models. The best 
studied one is the Nambu-Jona-Lasinio (NJL) model that was first developed to describe the 
interaction of nucleons |J and has later been applied for modeling low-energy QCD [p|-[TT 



with particular emphasis on the dynamical breaking of chiral symmetry and the occurence 
of the pion as a quasi Goldstone boson. The application of the NJL model for studies of 
quark matter thermodynamics is problematic since it has no confinement and free quarks 
appear well below the chiral phase transition [p~2] , [T3|1 . This contradicts to results from lattice 
gauge theory simulations of QCD thermodynamics where the critical temperatures for de- 
confinement and chiral restoration coincide. That can be helped by using a separable model 
which can be treated similarly to the NJL model but includes a momentum dependence for 
the interaction via formfactors. It has been shown [H| that in the chiral limit the model 
has no free quarks below the chiral transition. 

Looking again at the densities of compact star cores and comparing with the results of 
NJL model calculations |T^JT5[] it seems reasonable to include strange-flavor quarks in the 
model because the energy density is sufficiently high for their creation in weak processes. 
Therefore, we extend in the present work the separable model to the case of three quark 
flavors, assuming for simplicity U(3) symmetry. We will calculate the partition function 
using the method of bosonisation and applying the mean- field approximation. Finally we will 
formulate the resulting thermodynamics of three-flavor quark matter and obtain numerical 
results for the quark matter EoS and compact star structure. 



II. THE SEPARABLE QUARK MODEL 

The starting point of our approach is an effective chiral quark model action with a 
four-fermion interaction in the current- current form 
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r d 4 k 

+ J^k)j-(k')}] , (1) 

where we restrict us here to the scalar current jf(k) = q(k)X a f(k)q(k) and the pseudoscalar 
current j" (k) = q(k)ij 5 X a f (k)q(k) in Dirac space with q{k) and q(k) being quark spinors 
and the formfactor f(k) accounts for the nonlocality of the interaction. The action ([!]) 
is invariant under chiral rotations of the quark fields and color correlations are neglected 
(global color model). The generalization of previous models of this type |n, |TE| , |T7[ | to the 



three-flavor case is done by using the U(3) symmetry where A Q are the Gell-Mann matrices 
and A = y§J- 

Furthermore we do not include one of the possible models to account for the Ua(X) 
anomaly since, at least in the quark representation of the Di Vecchia - Veneziano model 
|T8| , it can be shown that there is no contribution to the quark thermodynamics on the 



mean-field level [jnj . For the quark mass matrix in flavor space we use the notation 

™ = E m / P /' ( 2 ) 
/ 

where rrif are the current quark masses and the projectors Pf on the flavor eigenstate 
f = u,d,s are defined as 

P u = ^A + \h + ^A 8 , (3) 

Pd = 7=6 X °-l X3 + ^ (4) 

Since the Matsubara frequencies in the T —>■ limit become quasicontinuous variables, 
the summation over the fourth component k± of the 4-momentum has been replaced by 
the corresponding integration. According to the Matsubara formalism the calculations are 
performed in Euclidean space rather than in Minkowski space where we use = The 
partition function in Feynman's path integral representation is given by 

= / VqDqexp \ S[q,q) - J ^^W 4 ^ , (6) 

where the constraint of baryon number conservation is realized by the diagonal matrix of 
chemical potentials fx (Lagrange multipliers) using the notation of the hat symbol analogous 
to©. 

In order to perform the functional integrations over the quark fields q and q we use the 



formalism of bosonisation (see |2(J and references therein) which is based on the Hubbard- 



Stratonovich transformation of the four-fermion interaction terms employing the identity 
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v n/^«pf^+/ t X*(*K(V) 



(7) 



for the scalar and a similar one for the pseudoscalar channel, where for the phase space 
integral the abbreviation J k = J has been used and M is a normalization factor. Now 
the generating functional is Gaussian in the quark fields and can be evaluated. We arrive 
at the transformed generating functional in terms of bosonic variables 

Z[T,fi]= JJVa a V7T a exp{S[a a ,7f a ]} (8) 

with the action functional 

S[a«, 7T«] = - J ^ In (det,^ + m + af(k) + ryte*/(*)]) 
a a a a n a 7i a 



with analogous use of the already known hat symbol and the 4- vector kf = [k,k4 + Z///J . In 
order to further evaluate the integral over the auxiliary bosonic fields a a and ir a we expanded 
them around their mean values a a and 7t a that minimize the action 

Oa = o a + o a {k) 

n a = TT a + 7l a (k) 

and neglect the fluctuations a a (k) and Ti a (k) in the following. The mean values of the pseu- 
doscalar field vanish for symmetry reasons ||19|| . The indices DFC refer to the determinant 
in Dirac-, flavor- and color- space. So we end up with the mean-field action 

S MF [T, {» f }) = E (-2N c j ^Pn(fcj + MJ)] + ^ , (10) 

with the effective quark masses My = Mf(k) = rrif + Aff(k) and the number of colors N c . 
The flavor dependent mass gaps Af are defined by a = J2f^fPf- 



A. Quark matter thermodynamics in mean field approximation 

In the mean field approximation, the grand canonical thermodynamical potential is given 

by 

n(T, {//}) = In {Z[T, {ii f }]/Z[0, {0}]} 

= {5 MF [T, {(if}) - 5 MF [0, {0}]} , (11) 

Where the divergent vacuum contribution has been subtracted. In what follows we consider 
the case T = only. In order to interpret our result, we want to represent it as a sum of 
three terms 
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n(o,{ji}) = 



(2vr 

+ y ( 2Nc [ ^-ln 
YV ■/ (2^) 4 



'A; 2 + (M°) 2, 



y y »- ^" J \ k 2 + m 2 

where M° = rrif + A®f(k 2 ) are the effective quark masses in the vacuum. The second 
term on the r.h.s. of this equation is the renormalized thermodynamical potential of an 
ideal fermion gas [^1]]. The third term of Eq. (|P2|) is independent of T and fi, i.e. it is a 
(thermodynamical) constant for the chosen model. Refering to the MIT bag model we call 
this term the bag-constant B. The remaining term includes the effects of quark interactions 
in the mean field approximation and can be evaluated numerically. 

All thermodynamical quantities can now be derived from Eq. (|12|). For instance, pressure, 
density, energy density and the chiral condensate are given by: 

„ n an _ on 

p V = -12 , n = -— , e = -p + fin , <q f q f >= ^— . (13) 

Still the quark mass gaps Af have to be determined. This is done by solving the gap 
equations which follow from the minimization conditions = 0. The gap equations read 

. r d 4 k 2Mff(k) , . 

A i = m{ -2N c )j mw Jin. (14) 

As can be seen from (|14|), for the chiral U(3) quark model the three gap equations for 
A u , Ad, A s are decoupled and can be solved separately. 



III. RESULTS FOR THE GAUSSIAN FORMFACTOR 



A. Parametrization of the model 



In the nonlocal separable quark model described above the formfactor of the interaction 
was not yet specified. In the following numerical investigations we will employ a simple 
Gaussian 



f(k) = exp(-k 2 /A 2 



(15) 



which has been used previously for the description of meson [22] and baryon [23|] properties 
in the vacuum as well as for those of deconfinement and mesons at finite temperature |I7| , [2"1 
A systematic extension to other choices of formfactors can be found in p5"tf26 



The Gaussian model has five free parameters to be defined: the coupling constant Do, the 
interaction range A, and the three current quark masses m u , m<j, m s . Setting m u = =: m q 
we restrict ourselves to four free parameters. These are fixed by the three well known 
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observables: pion mass m w = 140 MeV, kaon mass ttlk = 494 MeV and pion decay constant 
f n = 93 MeV. The formulas for the meson masses and the decay constant are calculated 
as approximations of the Bethe-Salpeter equation including the generalized Goldberger- 
Treiman relation [f26|1 . 



The fourth condition comes from values for the chiral condensate that are conform with 
phenomenology. The resulting parametrisations of the quark model are shown in Tab. |. 



B. Thermodynamics for quark matter without [3 equilibrium 

This case is relevant for systems which are considered for time scales larger than the 
typical strong interaction time of about 1 fm/c but smaller than the weak interaction time 
of several minutes, so that the presence of leptons (electrons) does not influence on the 
composition of quark matter and we can choose the chemical equilibrium with fi u = f^a = 
fi s = fi. For the numerical calculations we choose the parameter set for the light quark 
condensate —(uu + dd) 1 ^ 3 = 240 MeV which is a typical value known from phenomenology. 
We consider the behavior of thermodynamical quantities at T = with respect to the 
chemical potential. As we set m u = rrid earlier there is no difference between up- and 
down quarks and both are referred to as light quarks. Fig. [T] visualizes the behavior of the 
thermodynamical potential as a function of the light quark gap A q = A u = A^ for different 
values of the chemical potential fi. For fi < fi c = 333 MeV the argument and the value 
of the global minimum is independent of /i which corresponds to a vanishing quark density 
(confinement). At the critical value /i = fi c = 333 MeV a phase transition occurs from the 
massive, confining phase to a deconfining phase negligibly small mass gap. From the solution 
of the gap equation shown in Fig. ^ one can see that the strange quark gap of A s = 682 MeV 
still remains unchanged. Thus the strange quarks are confined until a higher value of the 
chemical potential fi C)S = 492 MeV is reached. This value is much bigger than the current 
strange quark mass. In Fig. |] we separate by vertical lines the regions of full confinement, 
two-flavor deconfinement and full deconfinement. Thus, in the present model, the onset 
of strange quark deconfinement is inhibited. Moreover, the onset of a finite strange quark 
density is not determined by a drop in the strange gap which remains constant and even 
starts to rise for large [i values. This result of the present model drastically differs from those 
of bag models or NJL models. The reason is the 4-momentum dependence of the dynamical 
quark mass function which results in complex mass poles for the quark propagators and 
makes the naive identification of the mass gap with a real mass pole impossible ||2S . 



The effect on thermodynamical quantities can be understood if we look at the pressure. 
In Fig. |3] we show for comparison the resulting equation of state for the pressure of the 
present separable model together with a two-flavor and a 3-flavor bag model. Both bag 
models are chosen such that the critical chemical potential for the deconfinement coincides 
with that of the separable model. The pressure of the present three-flavor separable model 
can be well described by a two-flavor bag model with a bag constant B = 81.3 MeV/fm 3 
in the region of chemical potentials 333 MeV < fi < 492 MeV where the third flavor is still 
confined. For comparison, the 3 flavor bag model has a bag constant B = 100.7 MeV/fm 3 
and is considerably harder than the separable one due to the additional relatively light 
strange quark flavor. 
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C. Inclusion of (3 equilibrium with electrons 



Quark matter in /3-equilibrium is to be supplemented with the two relations for con- 
servation of baryon charge and electric charge. In the deconfined phase there are quarks 
and leptons (in our model case up, down, strange quarks and electrons) with vanishing net 
electric charge 

2 1 

Q q (//\ //, //) + Q L (//) = -n u - -(n d + n s ) - n e = . (16) 

Taking into account the energy balance in weak interactions 

d <-> u + e~ + v e (17) 
s <-> u + e~ + v e (18) 

and introducing the average quark chemical potential // = + f^d + yU s ) we can write the 
/3 equilibrium conditions as 

2 1 
= /"--//, fi d = (i s = n + -n e . (19) 



Solving the equation of charge neutrality (16) one can find the chemical potential of electrons 
as a function of fi and using Eqs. ( |I9"D the equation of state can be given in terms of a single 
chemical potential fi. In Fig. f| we show the composition of the three-flavor quark matter 
for the Gaussian separable model in the case of f3 equilibrium with electrons as a function 
of the energy density. As for the case without (3 equilibrium we can define also in Fig. |] 
the regions of quark confinement (e < e c = 350 MeV/fm 3 ) and three-flavor deconfinement 
(e > e C)S = 930 MeV/fm 3 ). In the region of two-flavor deconfinement the concentrations of 
electrons, up- and down- quarks coincide with those of the two-flavor bag model except for 
the relatively small energies close to e c where the effect of a small dynamical quark mass 
leads to a density dependence of the composition. In Fig. § we demonstrate the influence 
of the P equilibrium on the equation of state. It can be seen that the difference between 
pressures with and without (3 equilibrium is limited to the region of intermediate densities, 
where the electron fraction reaches its maximum value x P ~ 0.002. 



IV. APPLICATIONS FOR COMPACT STARS 

One of the main goals for studying the strange quark matter equation of state is the 
possible application for compact stars. In particular, the hypothesis that strange quark 
matter might be more stable than ordinary nuclear matter |27j has lead to the investigation 
of possible consequences for properties of compact stars made thereof ]2"B|] . Most of these 
applications use the bag model equation of state where the result depends on the value 
of the bag constant as a free parameter. Recently, first steps have been made towards a 
description of strange quark matter within dynamical quark models such as the NJL model 
where the parameters are fixed from hadron properties. The non-confining quark 



dynamics of this model, however, leads to predictions for dynamical quark masses and 
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critical parameters of the chiral phase transition which differ from those of confining models 
3(|l4] and might be quantitatively incorrect. Here we want to extend previous studies 



of compact star properties with dynamically confining quark models to the strange quark 
sector and find the characteristics of stable compact star configurations with the equation 
of state derived above. 

For the calculation of the self-bound configuration for the quark matter with gravita- 
tional interaction one needs the condition of mechanical equilibrium of the thermodynamical 
pressure with the gravitational force. This condition is given by the Tolman-Oppenheimer- 
Volkoff- Equat ion 

dP , . . NN m(r) + 4irGr 3 P(r) , , 

- = -g W r) + P(r)) .I,,.;, _ 2gm(r) |' (20) 

and defines the profiles for all thermodynamical quantities in the case of nonrotating spher- 
ically symmetric distributed matter configurations in general relativity. In this equation m 
denotes the accumulated mass in the sphere with radius r given by 

FT 

m(r) = 47r / e(r')r' 2 dr (21) 



The gravitational constant is denoted by G. The radius R of the star is defined by the 
condition that the pressure becomes zero on the surface of the star P(R) = 0. The total 
mass of the star is M = m(R). 

Each configuration has one independent parameter which could be chosen to be e(0), 
the central energy density. In Fig. ^ we show the dependence of the total mass of the 
configuration as a function of the central density and the radius for the separable quark 
model and for the bag model in the cases of two and three flavors respectively. 

The rising branches of the mass-radius or mass-density relations correspond to the fami- 
lies of stable compact stars. The maximum possible mass for the separable model is 1.64 M Q 
for the three-flavor and 1.71 M & for the two-flavor case. The maximal central density is 
about 1350 MeV/fm 3 which allows for the three-flavor case to have strange quark matter 
in the core of the quark star. The comparison with the corresponding bag model strange 
stars shows that the latter are more compact, their maximum radius is about 8 km, and 
less massive with a maximum mass of about 1.5 M . The maximum radius of stars within 
the separable model is 11 km and thus exceeds the radii for both two and three flavor bag 
model quark stars. The origin of this difference is the behavior of the pressure in the low 
density region. 



V. CONCLUSION 

For neutron stars it is relevant to include the effects of strange flavor in a model for quark 
matter. In the simplest case considering £7(3) symmetry this can be done without increasing 
the complexity of the generating functional. We showed that in our separable model the 
gap equations decouple and can be solved separately. The resulting thermodynamics can 
be solved numerically and gives the equations of state for interacting quark matter. Unlike 
the well known NJL model the separable model is able to express the effects of confinement 
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in the thermodynamical quantities. The new result obtained within the present approach 
is the separation of the deconfmement of light quark flavors at fi c = 333 MeV from that 
of strange quarks which occurs only at a higher chemical potential of /x C)S = 492 MeV. A 
consequence for the application of the EoS presented here in compact star calculations is 
that strange quarks do occur only close to the maximum mass of 1.64 M Q , i.e. that for 
masses below 1.62 M , only two-flavor quark matter can occur. 
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TABLES 



^ — 7 7 „ 1 

[MeV] 


A 
[MeV] 


Do 
[GeV Z J 


m q 
[MeV] 


m s 
[MeV] 


a n 
[MeV] 


a n 
[MeV] 


230 


659.2 


29.32 


6.8 


143.5 


549.9 


767.8 


235 


697.6 


23.88 


6.4 


136.1 


497.0 


719.3 


240 


736.5 


19.88 


6.0 


129.5 


453.8 


682.1 


245 


775.5 


16.85 


5.6 


123.4 


419.7 


653.1 


250 


814.9 


14.49 


5.3 


118.0 


391.1 


630.0 


255 


853.8 


12.66 


5.0 


112.9 


368.7 


611.4 


260 


894.2 


11.11 


4.7 


108.1 


349.1 


596.1 



TABLE I. Parameter sets for the Gaussian separable model for different values of the chiral 
condensate < uu + dd > . 
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FIGURES 
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FIG. 1. Dependence of the thermodynamical potential on the light flavor gap A q 
(order parameter) for different values of the chemical potential, A s = 682 MeV. 
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FIG. 2. Solutions of the gap equations that minimize the potential 
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FIG. 3. Pressure of the quark matter as a function of the chemical potential for the separable 
model (solid line) compared to a three-flavor (dotted line) and a two-flavor (dashed line) bag model. 
All models have the same critical chemical potential fi c = 333 MeV for (light) quark deconfinement. 
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FIG. 4. Composition of three-flavor quark matter in [3 equilibrium with electrons. 
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FIG. 5. Pressure of three-flavor quark matter with (5 equilibrium and without. 
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FIG. 6. Stability for compact stars composed of quark matter. 
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